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The energy landscape of fullerene materials: a comparison between boron, 

boron- nit ride and carbon 
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Using the minima hopping global geometry optimization method on the density functional poten- 
tial energy surface we show that the energy landscape of boron clusters is glass like. Larger boron 
clusters have many structures which are lower in energy than the cages. This is in contrast to carbon 
and boron nitride systems which can be clearly identified as structure seekers. The differences in 
the potential energy landscape explain why carbon and boron nitride systems are found in nature 
whereas pure boron fullerenes have not been found. We thus present a methodology which can make 
predictions on the feasibility of the synthesis of new nano structures. 

PACS numbers: 36.40.Mr, 61. 46. Be 



The experimental synthesis of fullerenes is a very diffi- 
cult task. The carbon fullerene structures were therefore 
theoretically predicted [l[ long before they could be pro- 
duced in the lab [2[ . Many more hollow and enhodedrally 
doped fullerene structures made out of elements different 
from carbon have also been proposed since then theo- 
retically |3| in searches of other possible building blocks 
for nano-sciences. It is however surprising that since the 
experimental discovery of the carbon fullerenes some 25 
years ago no other fullerenes have been synthesized. So 
the question is whether experimentalists have just not 
yet found a way to synthesize these theoretically pre- 
dicted fullerenes, or whether they do not exist at all in 
nature. We have recently shown [4| that all the theoreti- 
cally proposed endohedral S120 fullerenes are meta-stable 
and can thus most likely not be found in nature. In this 
letter we investigate in detail boron clusters. Following 
the £?8o fullerene structure proposed by Szwacki et al. [5| 
various other fullerene [6[ and stuffed fullerene structures 
[?J were proposed. Subsequently it was however shown 
for Bgo that there exist non-fullerene structures [8| which 
are lower in energy. We will contrast the characteristics 
of the potential energy landscape (PES) of these boron 
clusters with those of systems found in nature, namely 
carbon and boron nitride fullerenes and find that there 
are important differences. 

To explore the energy landscape of the boron, carbon 
and boron nitride clusters we do global geometry opti- 
mizations on the density functional potential energy sur- 
face with the minima hopping algorithm [9j. This algo- 
rithm can render the global minimum configuration as 
well as many other low energy meta-stable structures. 
All the density functional calculations are done with the 
BigDFT electronic structure code [10j which uses a sys- 
tematic wavelet basis together with pseudopotentials [11] 
and the standard LDA Jll| and PBE [l2[ exchange cor- 
relation functionals. 

We start out by analyzing the BiqNiq cluster which 



was found to be short lived in experiments [l3[. In this 
system structural rigidity is imposed by a strong pref- 
erence for sp2 hybridization [14] as well as by the re- 
quirement that bonds are only formed between atoms of 
different type. This leads to a small configurational den- 
sity of states. As shown in Fig. [1] there exists a fairly 
large energy interval in which only cage like structures 
exist. Hence there is a strong driving force towards the 
ground state cage structure and minima hopping can find 
it rapidly. This driving force also allows the formation of 
£>16^16 in nature. 

Using the same methods we went on to study medium 
size boron clusters with 32-36 atoms. In this size range 
the clusters show a strong tendency to form cages and all 
the numerous low energy structures we found are cage 
like. This is agreement with a recent study [l5[ where 
the ground state was found to be cage like. Three rep- 
resentative ground state structures are shown in Fig. [2j 
These medium size clusters contain well known structural 
motifs [16| namely empty and filled hexagons as well as 
empty and filled pentagons. But in addition they contain 
numerous other structural motifs such as single atoms 
connecting filled hexagons or rings containing more than 
6 atoms. The inclusion of these other structural motifs 
does not rise the energy significantly and the first meta- 
stable structure is typically only 0.1 eV higher in energy 
than the global minimum. For .832 we found for instance 
some 100 cage like isomers in an energy interval of only 
1 eV above the global minimum and even more isomers 
presumably exist in this interval. The number of near- 
est neighbors in these structures varies from 4 to 6 and 
the bond angles vary from 90 degrees for some 4 fold 
coordinated corner atoms to 60 degrees for 6 fold coordi- 
nated atoms in the center of a planar hexagon. This is in 
contrast to the structural rigidity imposed by the sp2 hy- 
bridization on all the carbon fullerene structures we have 
generated. Even though one can find in our 60 atom car- 
bon structures rectangles and heptagons in addition to 
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FIG. 1: The configurational energy spectrum of BiqNiq. 
Boron atoms are shown in blue and nitrogen atoms in red. 
The higher energy cage structures can be described as a 'bas- 
ket' with a 'handle' made out of a chain of 4 atoms (two of 
each type). 
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FIG. 2: Global minima of ^32,^33 and B34 



hexagons and pentagons, all the atoms have, without any 
exception, 3 nearest neighbors in structures that are less 
than 20 eV above the ground state. As a consequence we 
expect the configurational density of states to be much 
smaller for carbon cages than for boron clusters. This is 
indeed the case as will be shown in more detail later. 
Next we did global geometry optimization runs for 



the £?8o cluster. A first run started from the Szwacki 
fullerene, which consists of the Cqq fullerene with 20 ad- 
ditional atoms filled into the hexagons. It thus consists 
of 20 filled hexagons and 12 empty pentagons. The inser- 
tion of the 20 atoms can be viewed as some kind of doping 
which stabilizes the two-dimensional boron network [17j . 
During a long period the cage structure was not destroyed 
in the minima hopping run. Instead minima hopping ex- 
plored the defect structures that we have described pre- 
viously [18| as well as other cage structures which are 
slightly lower in energy than the Szwacki fullerene. Since 
there is a very large number of possible defect structures 
this cage funnel contains a very large number of local 
minima and it takes long for minima hopping to escape 
from it. 

Once one escapes from the fullerene funnel one finds 
significantly lower energy structures. These structures 
contain the icosahedral B12 motif which is the basic 
building block of elemental boron. This icosahedron is 
in most cases at the base of a dome like structure or 
otherwise at the center of a spherical cage. Both the 
domes and the cages consist mainly but not exclusively 
of filled and empty hexagons and pentagons. Fig.[3]shows 
the configurational density of states for the B$q cluster. 
The majority of the structures are of the dome type and 
the energies of dome type and fullerene type structures 
overlaps. Like for the medium size boron clusters many 
structural building blocks can be combined to form clus- 
ters of very similar energy. Hence the energy difference 
between the low energy isomers is again very small. The 
lowest energy structure we found is considerably lower 
in energy than the recently proposed compact B 8 q struc- 
ture [8], both within the LDA and PBE functionals. 

Let us contrast the configurational energy spectrum of 
£?8o clusters with the one of Cqq clusters. For the first 
meta-stable structure is a Stone- Wales [l9| point defect 
which is 1.6 eV higher in energy than the fullerene ground 
state. Various defects can be combined to form cages of 
higher and higher energy. Two high energy structures 
are shown in Fig. 2J The lowest non-cage like structures 
are however some 25 eV higher in energy than the ground 
state. This shows that in contrast to Bgo the cage like and 
non-cage like structures are widely separated in energy. 
There is consequently a strong driving force towards cage 
like structures and finding the ground state for Cqo is 
much easier than for B$q. 

The differences in the potential energy landscape be- 
tween £?8o an d Ceo are also well illustrated by the follow- 
ing computer experiment. If one does a local geometry 
optimization for 80 boron atoms starting from random 
positions one obtains disordered structures which are al- 
ready fairly low in energy, namely about 10 eV higher 
than the ground state. This is in contrast to the case 
of 60 carbon atoms where a local geometry optimization 
starting from random positions gives structures which are 
about 50 eV above the ground state unless they happen 
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FIG. 3: The configurational energy spectrum of Bgo- The 
energy of the Szwacki fullerene is taken to be zero. The en- 
ergy levels of the icosahedron-dome structures are centered 
whereas the levels shifted to the left are fullerene like struc- 
tures. The levels on the right correspond to centered icosahe- 
dron structures. The atoms of the icosahedra are shown in 
yellow. The structure at an energy of -2.7 eV is the putative 
global minimum from ref [8|] . The energy per atom of our low- 
est energy B%o structure is about .13 eV per atom higher in 
energy than the sheet structure of Tang and Ismail-Beigi [I?} . 



to be cage like. This shows again that the boron poten- 
tial energy landscape has a glassy character with a lot of 
disordered low energy structures. The energy landscape 
of Cqo on the other hand has a broad and deep funnel 
which leads to the ground state fullerene. 

The glassy energy landscape of bulk boron has been ex- 
plained by the frustrated bonding features of boron where 
2-center bonds have to coexist with 3-center bonds [20|. 
The glassy energy landscape of the medium size boron 
clusters can also be explained in this way. Fig. [5] shows 
the coexistence of these two types of bonds in our lowest 
energy B 80 . 




FIG. 4: Two high energy Cqo cage structures. The structure 
on the left has only 3 fold coordinated atoms even though it 
contains two 7- member rings. It is 20.5 eV above the ground 
state. Structures that are even higher in energy can possess 
some chains with 2-fold coordination and anchor atoms for 
these chains with 4-fold coordination. The structure on the 
right is an example of such a cage and is 25 eV higher than 
the ground state. 




FIG. 5: The valence charge density in our lowest Bgo cluster 
(left) and the Cqo fullerene (right). Whereas in Cqo we see 
only two center bonds, both 2 and 3 center bonds are visible 
in B 8 o. 

In addition to the Bso cluster we also examined the 
B92 and £?ioo cluster. A structure with a icosahedron 
in the center of a 80 atom Szwacki fullerene is 7.7 eV 
lower than the fullerene which was obtained by filling 
the 12 pentagons [5]. The resulting structure has how- 
ever not anymore a high symmetry. A stuffed fullerene 
structure was proposed for £>ioo [2l| . Doing minima hop- 
ping runs starting from this configuration some struc- 
tures with lower energy and lower symmetry were found 
as well. These structures were also about 10 eV lower 
in energy than the recently proposed I?ioo fullerene [22| . 
This shows that disordered cages with an icosahedron in- 
side are the basic structural motif for boron clusters in 
this size range. 

Among all the ground state structures of boron clus- 
ters of any size, we could not find any high symmetries. 
Hence the vibrational modes have no or only low de- 
generacy. Following these modes by some mode follow- 
ing techniques will therefore in general lead to different 
transition states with different barrier heights. Since the 
height of the barrier correlates with curvature along the 
starting mode [23j , one can expect for a cluster of low 
symmetry a broader distribution of barrier heights and 
therefore a larger probability of finding low energy bar- 
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riers 24]. If low barriers exist a small modification of 
the external environment such as the presence of another 
cluster can make these barriers disappear. Hence it is 
not surprising that all boron structures that we exam- 
ined, independently of whether they are medium size, 
large, cage-like or not, turned out to be chemically reac- 
tive with other boron clusters when they are brought into 
contact. During such a chemical reaction with another 
cluster several chemical bonds are formed which leads 
to a considerable lowering of the energy and to a large 
distortion or even destruction of the original structures. 
This means that even though medium size clusters have 
a strong tendency for cage formation in isolation, it is 
unlikely that such boron cages exist in nature. This be- 
havior is also in contrast to the behavior of the Cqo and 
BigNiq fullerenes. They are only weakly interacting and 
do not form chemical bonds when they are brought into 
contact. The chemical reactivity of the boron clusters 
can also be rationalized in a local picture. If many differ- 
ent structural motifs can be used as a building block of 
a low symmetry cluster, it is very likely that some atoms 
have some dangling bonds which are chemically reactive. 

Our results explain why boron fullerenes have not been 
found experimentally. Boron clusters are frustrated sys- 
tems which do not have enough electrons to fill all elec- 
tronic orbitals in a chemical bonding based on pure sp2 
hybridization and they consequently do not exhibit some 
clear preference for a simple structural motif. Hence, 
from a energetical perspective, there is no driving force 
towards some well defined structure. Instead one finds 
a glassy energy landscape with a large number of differ- 
ent low energy structures whose energies are very similar. 
These structures are chemically reactive and will there- 
fore not be found under experimental conditions. The 
fact that no elemental boron but only compounds con- 
taining boron can be found on earth however indicates 
the possibility of synthesizing more complicated boron 
cages such as metal doped boron fullerenes. Such a dop- 
ing can energetically pull down the cage like part of the 
configurational space of boron clusters [18| . 

Our simulations demonstrate that one can make theo- 
retical predictions about the feasibility of an experimen- 
tal synthesis. In order to judge whether a system can be 
formed in nature , it is not necessary to simulate its syn- 
thesis process explicitly by molecular dynamics or sim- 
ilar methods. A global geometry optimization with the 
Minima Hopping algorithm indicates whether the system 
being simulated is a structure seeker or a system with a 
glass like potential energy surface. For a glassy system 
finding the global minimum is slow because one has to 
explore energetic regions with a large density of minima 
whose energies are very similar. For a structure seeker 
on the other hand the energy goes down rapidly and by 
significant amounts as one approaches the ground state. 
Only for these latter systems it is to be expected that 
synthesis pathways can be found. 



Our work thus clearly shows that theoretical cluster 
structure prediction has to be based on global geometry 
optimization because only this approach gives the neces- 
sary information on the potential energy landscape. The 
standard approach based on structures obtained from ed- 
ucated guesses, that were subsequently locally relaxed, 
gives only a very incomplete characterization of a sys- 
tem. A ground state structure predicted by global geom- 
etry optimization has a reasonable chance of being found 
in nature in significant quantities only if it is a) at the 
bottom of a broad and deep funnel, b) is significantly 
lower in energy than the other low energy meta-stable 
structures and c) has high symmetry. 
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